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ABSTRACT 

The Main Sequence Luminosity Function (LF) of the Galactic globular cluster NGC 288 has been 
obtained using deep WFPC2 photometry. We have employed a new method to correct for completeness 
and fully account for bin-to-bin migration due to blending and/or observational scatter. The effect of the 
presence of binary systems in the final LF is quantified and is found to be negligible. There is a strong 
indication of the mass segregation of unevolved single stars and clear signs of a depletion of low mass 
stars in NGC 288 with respect to other clusters. The results are in good agreement with the prediction 
of theoretical models of the dynamical evolution of NGC 288 that take into account the extreme orbital 
properties of this cluster. 

Subject headings: clusters: globular clusters: individual (NGC 288) - luminosity function 



1. INTRODUCTION 

After a decade of HST observations, the Luminosity 
Function (LF) of unevolved Main Sequence (MS) stars has 
been derived for several Galactic globular clusters (GGCs). 
In many cases the LF extends down to the faintest stel- 
lar objects (M ~ 0.1 Mq\ see, e.g., King et al. 1998; Be- 
din et al. 2001). For unevolved stars, the luminosity di- 
rectly tracks the mass, the most fundamental parameter 
of a star. The determination of the Initial Mass Function 
(IMF) from the observed LF of Main Sequence stars, as 
well as the study of the of the variation of the LF within 
a given cluster and among different clusters is now a well 
established field of research which has provided important 
insights on the dynamical evolution of GCCs (see, e.g., 
Ferraro et al. 1997; King et al. 1998; Marconi et al. 1998; 
Piotto & Zoccali 2000; Paresce & De Marchi 2000, and 
references therein). In particular, the signature of mass 
segregation, due to the energy equipartition established 
by two body encounters among clusters stars, has been 
detected in many clusters by comparing the LFs obtained 
at different distances from the cluster center. Also, the 
effect of the pruning of the outer extremities of the clus- 
ters by the Galactic tidal field or by the violent interaction 
with the Galactic bulge and disc have been shown with the 
analysis of the LF of MS stars (see, e.g., Andreuzzi et al. 
2001, and references therein). 

Here we present the MS LF of the globular cluster 
NGC 288 down to V ~ 24.5. The LF was based on 
HST-WFPC2 observations that are described in detail in 
a companion paper (Bellazzini et al. 2001b, hereafter Pap 
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I) , that was devoted to the study of the binary population 
of this cluster. In Pap I we used a method very similar to 
the one adopted by Rubenstein & Bailyn (1997), based on 
extensive sets of artificial star experiments, to estimate the 
global binary fraction (ft) of the cluster (7 %</(,< 30 % 
) and to show that binary systems are spatially segregated 
within the cluster. These same artificial star experiments 
can be used to obtain LFs in different regions of the clus- 
ter, corrected for all observational effects. In §2 we will 
briefly recall the characteristics of the adopted sample and 
of the artificial star experiments, as well as the results of 
Pap I that are relevant in the present context. In §3 we 
introduce a new method to correct the observed LF. The 
effect of the presence of binary systems on the final LF 
is also discussed and quantified with the support of the 
analysis performed in Pap I. In §4 we present the final LF 
of NGC 288 and compare it with the LFs of other clusters 
and, finally, in §5 we summarize our results. 

1.1. The globular cluster NGC 288 

NGC 288 is a loose ( logpo = 1-80 L QV pc~ 3 ; Djorgov- 
ski 1993) cluster. It is a classical old GGC (Rosenberg et 
al. 1999) with a blue Horizontal Branch (HB) morphology 
(see Bellazzini et al. 2001a; Catelan et al. 2001, and refer- 
ences therein) . Its stars are not exceedingly metal deficient 
([Fc/H] = —1.39) and have the abundance pattern typical 
of halo stars (Shetrone & Keane 2000). The cluster lies 
on a very inclined and eccentric orbit that carries it into 
the inner, denser region of the Galaxy, where the effects of 
Galactic tides and disc/bulge shocks are more disruptive 
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. According to Dinescu, Girard & van Altena (1999) the 
perigalactic distance is R p ~ 1.8 kpc, the eccentricity is 
e ~ 0.7, the orbital period is P ~ 220 x 10 6 years, and 
the "destruction rate" is one of the highest in the GGC 
system. The observed distance from the Galactic Center 
(Harris 1996, Rgc — 11 kpc;) indicates that the cluster is 
presently near the apogalactic point of its orbit (see also 
Pasquali, Brigas & De Marchi 2000). 

2. THE DATASET AND THE ARTIFICIAL STAR 
EXPERIMENT 

The observations, the data reduction process and the 
artificial star experiments are described in detail in Pap I. 
Here we summarize the points that are relevant in the 
present context. 

2.1. The sample 

Two HST-WFPC2 fields were observed in the F555W 
and F814W passbands: one with the PC camera nearly 
centered on the center of the cluster (Int field), and one 
with the WF4 camera partially overlapping the WF3-Int 
field (Ext field; see Fig. 1 of Pap I). The Int field is almost 
completely enclosed in the region within 1 half light radius 
(1 rh — 1-6 r c ; where r c = 85 is the core radius) Djorgov- 
ski 1993) while the Ext field samples a region contained 
in the annulus lr/, < r < 2 r^ , where r is the distance 
from the cluster center. The data were transformed to 
the standard Johnson-Kron-Cousins photometric system 
by calibration to the ground-based data of Rosenberg et al. 
(2000). The Int sample contains 5766 sources from V ~ 13 
to V ~ 25 and the Ext one contains 2013 sources between 
V ~ 16 and V ~ 25.5. The large majority of the observed 
stars are fainter than the Turn Off (TO) point (V ~ 19), 
and hence are unevolved MS stars (and/or binary systems 
containing two unevolved MS stars; see Pap I). The F336W 
and F225W observations of the Int field that has been 
shortly discussed in Pap I are not used in the present pa- 
per. 

2.2. The artificial star experiments 

Simulations containing more than 80,000 artificial stars 
per chip have been made on the V, I frames. In each run 
of the artificial stars experiments ~ 100 artificial stars per 
chip were added at the same position in the median image 
(that we used as a master-frame to obtain a list of posi- 
tions of bona fide stars; see Pap I) and in each of the single 
frames (that were used to obtain accurate photometry; see 
Pap I) . The set of experiments for a given chip (sub-field) 
was completed with ~ 800 independent runs. The final 
total number of artificial star experiments is > 1, 500, 000. 
For each artificial star a random input V magnitude is 
drawn from a distribution similar to the observed LF ex- 
trapolated to magnitudes slightly below the observed lower 
limit. This fainter limit was required so that the simula- 
tions would contain a large number of stars fainter than 
the observed limiting magnitude to correctly sample the 
magnitude bins where the incompleteness is expected to 
be highest (see Pap I). 

We ultimately derive the true LF in a way that does not 
depend on the assumed artificial star LF (see §3.2). The 
input I is obtained from the input V using the cluster 
ridge line. Each simulated star thus had the appropriate 



color and was simultaneously added to the corresponding 
V and I frames at the same position. The simulated stars 
were randomly distributed on the frames with an addi- 
tional constraint (see Pap I and also Tosi et al. 2001, for 
a detailed description of the method) preventing any in- 
terference between simulated stars. In practice, for each 
run of the artificial stars experiment, only one artificial 
star is added in each 80 x 80 pixel portion of a chip. The 
whole process of data reduction has been repeated on all 
the frames with artificial stars added and the output mag- 
nitudes have been recorded (see Pap I, for details). 

Fig. 1 shows the input and output CMDs for the stars 
simulated in the PC-Int, PC-Ext, WF2-int and WF2-Ext 
frames. The WF2 samples have been chosen as repre- 
sentative of the WF cameras for the Int and Ext fields. 
Within a given field (Int and Ext) the observations with 
the different WF cameras have very similar properties (see 
Pap I). Note the different response to the same input color- 
magnitude distribution in the four cases represented. This 
arises because the larger pixels of the WF cameras allow 
them to reach fainter limiting magnitudes than the PC 
camera and because the Int and Ext fields have different 
average degrees of crowding. 

All of the effects of the whole process of observation 
and data reduction on a large sample of stars represent- 
ing a pure Simple Stellar Population (SSP) are shown in 
Fig. 1. These effects can be imagined to be the results of 
the application of a multi-dimensional function (that we 
will call Observation/Measure Function; hereafter OMF) 
on the true stars of the clusters, here represented by the 
input artificial stars lying on the ridge-line. In Pap I the 
empirical knowledge of the effects of the OMF has been 
used to extract the signal provided by the presence of a 
significant population of binary systems in the color distri- 
bution around the ridge-line from the observational noise 
due to random scatter, blending, etc. Here we will use it to 
correct the observed LF from all the spurious deformations 
due to the same factors. 



2.3. Completeness and division in subsamples 

The artificial star experiments have been carried out in- 
dependently for each camera of each observed field. The 
completeness functions for each camera have been pre- 
sented in Fig. 5 of Pap I. The crowding conditions of the 
dataset are never particularly critical, and the complete- 
ness is rather similar everywhere. The completeness factor 
(Cf) is larger than 80 % for V < 23.5 in all the samples 
except the PC-Int, the innermost one, where Cf > 80 % 
for V < 22. It has been shown that the spatial variation 
of completeness within a single camera are negligible. For 
the present purpose, we take advantage of the great sim- 
ilarity of crowding conditions in the WF cameras of the 
Int and Ext samples to subdivide the total sample in the 
following subsamples: (1) PC-Int, (2) WF-Int, contain- 
ing the WF2-Int, WF3-Int and WF4-Int samples, (3) PC- 
Ext and, (4) WF-Ext, containing the WF2-Ext, WF3-Ext 
and WF4-Ext samples. The four subsamples are internally 
very homogeneous and cover different radial regions of the 
cluster. In this sense, the difference between the PC-Ext 
and the WF-Ext samples is small but we prefered to keep 
them separate because of the different limiting magnitude. 
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Fig. 1. — CMDs of the simulated stars obtained with the input magnitudes (Vi n , 7i n ; left panels), and with the output magnitudes (V ou t, /out! 
right panels). The plot shows the effect of the observation - reduction process on the PC-Int and PC-Ext frames, while WF2-Int and WF2-Ext 
are shown as representative of the WF cameras in the Int and Ext samples, respectively. Ns is the total number of simulated stars while 
Nr is the total number of recovered stars. Note that the ratio between this two numbers is not directly comparable, at face value, since 
the limiting magnitude of the LF from which the artificial stars were drawn was different in the different cases. Since the observed limiting 
magnitude was different in the various cases, the limiting magnitude of the artficial stars LF was varied accordingly to exclude unnecessarily 
faint stars (e.g. more than ~ 1 mag below than the observed limiting magnitude) from the simulations. 
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3. THE CORRECTION OF THE OBSERVED LF 

The effects of the OMF on the true properties of stars 
that may affect the LF derived from an observed sample 
are: 

1. Star loss. Stars may be lost in the observa- 
tion/measure process. This occurs for two different 
reasons: 

(a) Near the limiting magnitude of the observa- 
tions a star can be succesfully detected if it 
falls on a positive fluctuation of the back- 
ground that brings it above the detection 
threshold and can be lost if it falls on a nega- 
tive fluctuation. This effect is purely stochas- 
tic; it does not depend on the crowding. It just 
depends on the amplitude of the random fluc- 
tuations of the background. There is a level 
of magnitude beyond which no possible posi- 
tive fluctuation can bring a star above the de- 
tection threshold. All the stars fainter than 
this limit are lost (unobserved) and just con- 
tribute to the diffuse luminosity of the unre- 
solved background. 

(b) Stars may be lost if they fall on a brighter 
source (star or galaxy) and/or if their images 
are badly corrupted by defects of the chip or 
the hit of a cosmic ray (this two latter effects 
are expected to have a negligible impact in the 
present case because of the adopted observa- 
tion + data reduction strategy; see Pap I). 
This effect depends on the local degree of 
crowding as well as on the magnitude of the 
star. 

2. Blending. A star may be blended with a fainter 
source and its measured magnitude will be brighter 
than the true one because of the flux contributed 
by the fainter blended source. Considered from the 
perspective of artificial stars experiments the max- 
imum possible effect of a blending is to change the 
observed magnitude of a star by ~ 0.75 mag. In 
fact, if we recover an artificial star by its (known) 
position in a frame with an output magnitude more 
than 0.75 mag brighter than its input magnitude, 
this means that it is fallen on a real star brighter 
than the simulated star. Thus, we are not recover- 
ing/measuring the simulated star but the real one 
instead, i.e., thus the artificial star has been lost. 

3. Observational scatter. The magnitude of the 
true star is altered by all the sources of stocas- 
tic noise associated with the process of observa- 
tion/data reduction (e.g., photon noise, imperfect 
PSF fit, etc.). The maximum possible effect of this 
factor is magnitude dependent and may be roughly 
evaluated on the basis of the typical error on a given 

4 The error associated to each magnitude entry is the error on the mean 

5 Note that this phenomenon may significantly alter the star counts in 
crowding conditions are critical, see Tosi et al. (2001), for discussion 

6 The input magnitude (VS n ) of artificial stars corresponds to the true 
corresponds to the observed one (V j,,,) under the extant conditions 



measure. For instance, in the present case we in- 
cluded in the final samples only the stars with er- 
rors 4 lower than 0.2 mag in both / and V. Hence 
a ~ 3<t fluctuation can change the magnitude of an 
observed star by ~ 0.6 mag, at most, near the lim- 
iting magnitude level. 

The derived LFs are in practice histograms of the num- 
ber of stars as a function of magnitude. The factors above 
change the true LF histogram by depopulating the bins 
as a function of (increasing) magnitude (mainly factor 1) 
and by moving stars from one bin to another (migration) 
due to the magnitude changes induced by factor 2 and/or 
3. However factors 2 and 3 are strictly interwoven and 
their effect cannot be separated. For example, a given 
star may be moved, say, 0.4 mag brighter because it is 
blended with a fainter unresolved source, but it might si- 
multaneously be affected by a negative fluctuation of the 
random observational scatter, so its final magnitude may 
differ from the true one by just 0.2 mag. A blending be- 
tween two unresolved sources can bring them above the 
detection threshold, injecting in the sample two stars that 
should have been lost if not blended 5 . The possibilities 
are too numerous to list here. The real contribution of the 
three effects cannot be known for a the specific case. On 
the other hand, it is clear that we can obtain a statistical 
knowledge of the actual LF. 

The traditional way to recover the true LF from an ob- 
served one consists in multiplying the numbers in each bin 
of the observed LF by the inverse of the Completeness fac- 
tor (Cf = Nr/Ns; where Nr is the number of recovered 
stars and Ns is the number of the simulated ones) at the 
center of the bins. It is clear that such procedure corrects 
at first order only for the effects of star loss (factor 1) 
while the effects of the other factors are ignored. A popu- 
lar misconception is that the observational scatter cannot 
produce bin migration if the width of the bin is larger than 
the typical error on the estimate of the magnitude. This 
is obviously not true since a star whose true magnitude 
is near the edge of a bin can migrate because of an arbi- 
trarily small fluctuation. The effects of migration may be 
amplified by the completeness correction, especially in the 
faintest bins of the LF. 

Furthermore the Completeness factor may be obtained 
from the artificial star experiments only as a function of 
input magnitude (say, Vi n ) while the observed LF is (obvi- 
ously) a function of the observed magnitude, whose corre- 
sponding quantity in the artificial star experiment is V^ut 6 . 
Thus when an observed LF is corrected by multiplying it 
by l/C/(Vin), it is implicitly assumed that Vi„ and V b s 
are the same quantity, that, in general, is not true. 

There have been few attempts to correct the observed 
LFs for all the effects above (Drukier et al. 1988; Stetson 
& Harris 1988). The proposed methods require a complex 
mathematical approach, and may rely on the a priori as- 
sumption of an analytical form of the true LF (Stetson & 
Harris 1988) or on approximate treatments of the propa- 
gation of errors under matrix inversions (see Drukier et al. 

over repeated measures, see Pap I. 

the bins near the limiting magnitude, at least in cases in which the 
magnitude of real stars (Vtrue), while the output magnitude (Vout) 
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1988; Mighell 1990, for references and discussion). 

Here we propose a new method that is fully non- 
parametric and extremely simple in nature, taking also ad- 
vantage of the huge enhancement in computational power 
of workstaions and PCs that took place since the times of 
the above quoted attempts. As we will show below (a) the 
proposed method takes full account of the effects of bin mi- 
gration and (6) the completeness correction is consistently 
applied to a LF that is function of V\ n as the completeness 
factor. The method is particularly appropriate for the 
application to LFs of SSPs, but it can be probably gener- 
alized to other cases. The prerequisite is the availability of 
huge sets of realistic artificial star experiments well sam- 
pling the whole spatial, luminosity, and color range of the 
observations. 

3.1. The Equivalent Sample 

Consider a large set of artificial star experiments per- 
formed in a region of uniform crowding conditions, e.g., 
the WF2-Ext set presented in Fig. 1. Suppose that the 
input magnitudes of the artificial stars be distributed as 
the "true" stars of the parent population (in this case the 
stars of NGC 288; the reasons of this hypothesis will be- 
come clear at the end ofthis section). From this set it is 
possible to extract a subsample of recovered stars with ex- 
actely the same dimension and the same distribution as 
a function of V ou t as the observed sample. Call this sub- 
sample an Equivalent Sample (ES) 7 . It is important to note 
that the ES contains only stars that have been succesfully 
recovered in the artificial star experiments process, hence 
that have suffered only the effects of bin migration (i.e, 
blending and/or observational scatter). 

By definition the ES has the same output ( e.g., observed) 
properties of the observed sample, it is one of the many 
possible relizations that are mapped by the OMF from 
the input/true domain to the output/observed one. Thus 
the input (e.g., true) properties of the ES are one of the 
possible true set of properties that may have been mapped 
into the observed sample by the OMF. It is immediately 
apparent that the LF of the input magnitudes of the ES 
is equivalent to the observed LF once corrected for all the 
effects of bin migration® . Thus the identification of an ES 
allows (a) the removal of all the effects of bin migration 
and (b) the expression of the observed LF as a function of 
input magnitude, consistently with the completeness fac- 
tor (sec §3.). If we apply the completeness correction to 
the input LF of the ES we find one possible realization 
of the final LF, corrected for all the observational effects. 
Since different ES may be mapped into the same observed 
sample by the OMF, repeating the process for a number 
of (randomly extracted) different ES and taking the aver- 
age will allow to obtain a final LF more robust to random 
fluctuations as well as a direct estimate of the uncertainty 
associated with the whole process. 

However, there is another important factor we have to 
take into account. While the completeness correction is 
completely independent from the distribution in magni- 
tude of the artificial star set, the derived correction for 



bin migration is not. Consider a given bin of the ob- 
served LF. Such a bin will be primarily filled by stars with 
y out ~ Vi n , but also by a certain fraction of stars which 
have migrated from the nearby bins beacause of blending 
and photometric errors. Increasing the number of artifi- 
cial stars in the nearby bins, with respect to the considered 
bin, will enhance the probability of extracting stars that 
enter in the bin because of migration, hence increasing the 
final "migration correction" . The right correction would 
be obtained only if the artificial star set is distributed as 
the true LF of the population which is unknown, being the 
very target of our analysis. We circumvent this problem 
by iteratively adjusting (a posteriori) the distribution of 
the artificial star set until the best approximation is ob- 
tained (see below for details). The approach has proven to 
be successful since we obtain the same final LF indepen- 
dently of the initial distribution of the artificial star set. In 
the following section the adopted operational procedure is 
described in detail. 

3.2. Operational procedure 

For each of the four observed subsamples (PC-Int, WF- 
Int, PC-Ext, WF-Ext) we adopted the following proce- 
dure: 

1. An initial distribution of input magnitudes is as- 
sumed for the artificial stars and the largest possible 
subset accordingly distributed is extracted from the 
whole set of artificial stars. 

2. Twenty different ES are extracted at random from 
the above subset and their LF(Vi n ) are corrected for 
completeness. 

3. The (bin per bin) average of the twenty obtained 
LFs is adopted as the "current corrected LF" and 
the standard deviation is adopted as the correspond- 
ing uncertainty. 

4. The "current corrected LF" is assumed as the dis- 
tribution of the artificial stars and a new subset of 
artificial stars accordingly distributed is extracted 
from the whole set of artificial stars as in step 1. 

5. The steps 2,3 and 4 arc repeated twenty times. The 
twentieth "current corrected LF" is adopted as the 
final LF with the corresponding uncertainties. 



Step 4 ensures that at each new iteration the adopted 
set of artificial stars will become more and more similar 
to the "true" one, giving increasingly appropriate rates of 
bin migration. 

In all the cases considered here the process converged in 
< 6 iterations, e.g., the "current corrected LFs" become 
very stable after the sixth iteration. Furthermore, the re- 
sult is independent of the initial distribution of input mag- 
nitudes that is assumed for the artificial star subset. We 
have tried many different distributions (some of which are 

7 The operational procedure to extract an ES is the following. Consider the histogram representing the observed LF: each bin i contains N{ 
stars. From the artificial star set we extract at random Ni stars whose output magnitudes lie in the i-th bin. Repeating this step for each 
bin of the observed LF we obtain a sample of recovered artificial star that have an output LF equal to the observed one. Note that the same 
procedure can be applied on a star-by-star basis, by randomly extracting an artificial sister (see §3.4) for each observed star. 

8 In other words the LF of the input magnitudes of the ES is the observed LF seen before any blending and/or observational error take place. 
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shown in Fig. 2) and in all cases the process converged to 
the same final LF in few iterations. 

Some examples of the procedure are shown in Fig. 2, 
applyied to the WF-Ext sample. In Fig. 2a the artificial 
star subset was initially distributed as the observed LF 
down to V < 25.5 and uniformly from this limit down to 
V = 27. A scaled version of the adopted distribution is 
plotted as a thick dotted line. The thin dotted histogram 
shows the "current corrected LF" after the first iteration 
of the process, the thin continuous histogram shows the 
"current corrected LF" after the sixth iteration and the 
open dots are the final LF, corresponding to the twentieth 
iteration of the process, with the currrent standard devi- 
ation of the twenty realizations of the corrected LFs as 
errorbars. The encircled dot is the first bin of the LF for 
which Cf < 0.5, i.e. the completeness correction is larger 
than a factor 2. The fast convergence and the stability of 
the final solution can be readily appreciated. 

Fig. 2b and 2c have the same symbols of Fig. 2a but 
shows the convergence of the process starting from (very) 
different initial distributions of the artificial star subset. 
The open dots are the same final LF shown in Fig. 2a. 
Thus Fig. 2a and 2b shows that the process converges to 
the same final LF independently of the assumed initial 
distributions of the artificial stars subset. 

Finally, Fig. 2d shows in detail the results of the twen- 
tieth (i.e. final) iteration of the process shown in Fig. 2a. 
The twenty histograms are the corrected LFs obtained by 
the twenty different ES extracted from the finally adopted 
subset of artificial stars. The full dots and errorbars are 
their means and standard deviations. 

3.2.1. The final LFs 

Since the obtained LFs are independent from the ini- 
tial distribution of the adopted subset of artificial stars, 
we derived our final (corrected) LFs for all the considered 
samples assuming the same initial distribution of the ar- 
tificial stars shown in Fig. 2a. The final LFs are given 
in Table 1 (PC-Int and WF-Int samples) and in Table 2 
(PC-Ext and WF-Ext samples). The table gives the V in 
(e.g., Vtrue) of the center of the bins and, for each sample, 
the observed LF, the completeness factor C/, the final LF 
and the total error on the final LF (see below). 

The final LF is given down to the first bin that has Cf < 
0.50. In the following plots this point is circled to serve as 
a reminder that it is has been derived applying a large cor- 
rection. We retain these points as educated guesses of the 
behaviour of the LFs beyond the range in which moderate 
and safe corrections can be made. We consider as fully 
reliable only the points of the LF for which Cf > 0.50. 
We note that the last fully reliable points, according to 
the above criterium, have C f = 0.607, 0.769, 0.806, 0.842 
for the PC-Int, WF-Int, PC-Ext, and WF-Ext samples, 
respectively. Thus the completeness is quite high for the 
whole considered range. 

The final total error reported in Tab. 1 and 2 has been 
obtained by summing in quadrature the error on the Com- 
pleteness factor (estimated according to Eq. 18 by Mighcll 
1990), to the standard deviation of the twenty different 
realizations of the LF that are obtained in the last step of 
the procedure (see §3.2 and Fig. 2d). In all the following 
plots the errorbars associated with the LFs of NGC 288 



represent this final total error. 

The final LFs (open dots) are shown with the corre- 
sponding observed LF (full dots) in Fig. 3. 

As a further check we compared our final LFs with those 
obtained by dividing the observed LF(V b 3 ) by C/(Vi n ). 
We note that the differences are within lc in most bins 
and < 2cr in any bin. As expected, in this particular ap- 
plication the effects of bin migration are not large. This 
would not be the case if much more crowded fields were 
considered (see §3. and Tosi et al. 2001). 

3.3. The effect of binary systems on the LF 

There is an intrinsic kind of blending that cannot be 
directly tackled with artificial star experiments. Physical 
binary systems are strictly equivalent to chance superposi- 
tion blcndings. The LF for single stars cannot be obtained 
if an independent estimate of the incidence of binary sys- 
tems is not available. Further, since the actual binary frac- 
tion (/),) depends on the (unknown) distribution of mass 
ratios (F(q); see Pap I), any correction would be model 
dependent. Indeed, a direct estimate of the impact of bi- 
nary systems on the LF of a globular cluster has only been 
possible for NGC 6752 (Rubcnstein & Bailyn 1999). 

Now we can add a second cluster using the results of 
Pap I. We measured fb in the Int and Ext regions of 
NGC 288. In the Int region we found that the observa- 
tions are compatible with 8%< fb < 38%, and the most 
probable range is 10% < fb < 20%. In the Ext region < 
10 % and the most probable value of fb is zero, indepen- 
dently of the assumed F(q). Thus, a sizeable fraction of 
binary systems is present in NGC 288, at least in the Int 
sample, and may affect the derived LF of single stars. 

To quantify the effect we computed the mean bin-by- 
bin ratio between the LFs obtained by 5 ES with fb = 20 
% and by 5 ES with fb = %. This number gives (ap- 
proximately) the fractional variation of the star content of 
the bins due to binary systems. Furthermore, dividing the 
final LF by the quoted ratio, the LF is corrected for the 
binary content and mapped to the "pure single stars" case. 
In Fig. 4 the final LF of the WF-Int sample not corrected 
for the binary content (filled dots with error bars) is com- 
pared with its binary-corrected versions in the assumption 
fb = 20 % and for three different F(q) (see caption and 
Pap I for details). The comparison is performed in the 
magnitude range in which the binary estimate of Pap I 
has been performed (20 < V < 23.5). 

There are several noteworthy features in Fig. 4: 1) The 
maximum correction is < 3%. 2) All the corrected LFs 
lie within the uncertainties of the un-corrected one. 3) 
The PLMR F(q) case is most similar to the uncorrected 
case. This was expected, since in this case the large ma- 
jority of binary system has a very low mass secondary. 
This means that the migration of the primary is small, 
and that most of the secondary stars are so faint that 
they fall in an unobserved region of the LF. 4) The cor- 
rection for the PHMR case depletes the intermediate bins 
(21.5 < V < 23; because it corrects for the migration of 
primary) and enhances the last bin because of the recov- 
ery of some secondary stars as single stars. 5) No bin is 
exclusively depleted or enhanced by the redistribution of 
the binary components as single stars. 

From the above results and discussion we conclude that: 
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Fig. 2. — Panel (a): The final LF is shown as open dots with error bars (the same LF is reported with the same symbols also in panels b 
and c). The encircled dot is the first point with Cj < 50 %. The thick dotted line is a scaled version of the distribution of Vi n of the artificial 
stars subsamplc adopted for the first iteration, in the present case assumed with the same shape of the observed LF to V = 25 and with 
N = 215 for the V > 25 bins (in the present scale). The thin dotted histogram is the LF recovered after the first iteration of the process. The 
thin continuous histogram is the LF recovered after the sixth iteration, when convergence has been reached. Panel (b): the same as Panel (a) 
but adopting a different distribution for the artificial stars subsample adopted for the first iteration (i.e., flat for V > 20.5). Panel (c): the 
same as panels (a) and (b) for an initial distribution for the artificial stars subsample strongly peaked at V = 22.75. Panel (d): the twenty 
realizations of the last iteration of the process for the case shown in Panel (a). The final LF is the mean of the twenty realizations (full dots) 
and the uncertainty associated with the process is the corresponding standard deviation. 
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1. The correction for binary systems can be safely ne- 
glected for the Ext samples where the binary frac- 
tion is much lower than f b = 20 % and possibly 
null. 

2. The correction can be neglected also for the Int sam- 
ples, at least for V < 23.5, since it is smaller than 
the uncertainties of the final LF. However, it cannot 
be excluded that the faintest bins of the single stars 
LF need a significant enhancement to take into ac- 
count the effect of binaries, particularly if the actual 
F(q) favors high mass ratios. 

3.4. Further developements and other applications 

The case under study is particularly simple. NGC 288 is 
a very loose and relatively nearby cluster, so our WFPC2 
observations well resolve the large majority of the stars 
down to the limiting magnitude of each subsample with 
a very high degree of completeness. Further, the density 
gradient is a very fortuitous match to WFPC2 allowing a 
few large subsamples to nicely probe radial variations. 

It is important to note that the Equivalent Sample con- 
cept does not require such simplifications and may be suc- 
cessfully applied under any condition. If a strictly rigorous 
approach is adopted the ES approach may have far reach- 
ing applications in the study of stellar populations. Here 
we provide two possible examples to show the potential of 
the method. 

• Extremely crowded field with strong density 
gradient. The method can tackle with samples 
with arbitrarily large density gradient, for instance 
by subdividing the observed sample in subsamples 
whith nearly uniform stellar density (as done in the 
present application). If the dimension of the ob- 
served sample is not sufficiently large to allow this 
approach, the ES method can be implemented also 
on a star-by-star basis. It will suffice to assemble 
the ESs by associating an artificial star analogue for 
each observed star, randomly extracted from a set of 
artificial stars having (approximately) the same V^ u t 
and position in the frame of the considered observed 
star (artificial sisters). The only basic requirement 
to successfully follow this approach is to have a very 
large set of artificial stars, sampling with large mul- 
tiplicity the whole range of magnitudes and posi- 
tions covered by the observed sample. 

• Mapping synthetic CMDs into the observa- 
tional plane. The standard tool used to derive the 
Star Formation History (SFH) from the CMD of a 
resolved galaxy consist in the reproduction of ob- 
served CMD by a synthethic CMD in which the stars 
are extracted from theoretical evolutionary tracks 
(see Lejeune et al. 2001, for references and discus- 
sion of the various adopted techniques). To trans- 
form the synthetic CMD to the observational plane 
the effects of incompleteness, blending and observa- 
tional scatter must be added to the synthetic stars. 
This task is not easy and it is usually done with sim- 
plified approaches. For instance, the observational 
scatter is usually introduced by randomly extracting 
'errors' from gaussian distributions with standard 



deviations equal to the typical la photometric er- 
ror. However, the distribution of the observational 
scatter is not gaussian, in general, and it is not easy 
to parametrize. Furthermore, the effects of blending 
are not separable from the effects of observational 
scatter, and gaussian distributions do not provide 
a good reproduction of the real distribution, intro- 
ducing undesircd noise in the comparisons. 

With the ES approach, one can easily map the input 
magnitudes (e.g., "true" magnitudes) extracted by 
the evolutionary track in to the observational plane 
with the correct OMF, by "asking" to the artificial 
stars experiments what is the probabilty that such 
star would be successfully recovered or not. If a 
star passes the "incompleteness barrier", it can be 
correctely mapped into the observational plane, for 
instance by assigning to it the output properties of 
an artificial analog extracted at random from a set 
of artificial sisters. 

In this way, the whole effects of the OMF would be 
correctly reproduced without dangerous approxima- 
tions and/or parametrizations. 

4. COMPARISON OF LFS 

Once the observed LFs have been corrected for all the 
observational effects and the impact of binary systems 
have been quantified, we can use the final LF to study the 
dynamical status of NGC 288 either comparing the LFs 
in the Int and Ext regions of the cluster, either comparing 
the LF of NGC 288 with those of other clusters. In the 
following we will normalize all the samples to be compared 
to the number of stars in the range 18.5 < V < 20.5, fol- 
lowing the approach of Piotto & Zoccali (1999, hereafter 
PZ99). 

NGC 288 is not sufficiently nearby to allow a safe deriva- 
tion of the LF down to the hydrogen burning limit even 
with HST. The comparison with other "state of the art" 
LFs of more nearby clusters (see Fig. 6, below) shows that 
such limit would occur a couple of magnitudes below our 
faintest valid point (the faintest bin of the WF-Ext LF) . 
The Mass-Luminosity relation adopted in Pap I indicates 
that the lowest mass efficiently sampled by our deepest LF 
is M ~ 0.35 Mq. Usually MF are compared according to 
their slope in the plane logiV — logM for M < 0.5 M Q 
(see PZ99 and references therein). Our best LF has just 
three points below this limit, thus the slope of the MF 
would be poorly constrained. Because of this limitation 
we avoid any conversion of our LFs to Mass Functions 
(MF), since that would necessarily be somewhat model 
dependent and uncertain (see, e.g., Bcdin et al. 2001). In- 
stead we will make direct comparisons between LFs in the 
allowed range. 

4.1. Mass segregation within single stars in NGC 288 

The LFs of the PC fields are identical (to within the er- 
rors) to those of the respective WF samples, in the range 
where the comparison is possible. Since the WF samples 
are much larger and reach fainter magnitudes we will limit 
our analysis to these homogeneous datasets in the follow- 
ing. 

In Fig. 5 the LFs from the WF-Int and the WF-Ext sam- 
ples are compared. It is evident at a first glance that the 
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Fig. 5. — Comparison between the LFs of the WF-Int sample (open circles) and of the WF-Ext sample (open triangles). The encircled 
points are the first points with Ct < 0.50. The horizontal bars indicate the normalization interval. 
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LF of the WF-Int sample is flatter than that of WF-Ext. 
Given the adopted normalization, this is a clear indication 
of a depletion of low mass stars in the region within 1 
(Int) with respect to the Ext field. If the brightest bin 
is included in the normalization the effect is significantly 
enhanced. We note also that any correction for binaries 
would be much smaller than the observed difference (see 
Fig. 4). On the other hand, if fainter stars (20.5 < V < 24) 
are used for normalization, one would infer a substantial 
excess of the more massive stars in the Int sample. 

Independently of the assumed normalization it is worth 
establishing if the difference between the two LFs is statis- 
tically significant. To check this hypothesis we performed 
a x 2 test, adopting the appropriate definition of the statis- 
tic according to Eq. 14.3.3 by Press et al. (1992). The 
comparison has been performed in the safely comparable 
range 18.5 < V < 24. The resulting reduced x 2 is 2.31 
with 12 degrees of freedom. The difference between the 
two LFs is found to be significant at the 99.1 % confidence 
level, i.e. highly significant. 

We have already demonstrated that mass segregation 
exists in NGC 288, since both binary systems and their 
by-product Blue Stragglers Stars are also more centrally 
concentrated than the the single stars. Although the evi- 
dence is not as strong (because of the possible ambiguity 
in the interpretation associated with the choice of the nor- 
malization range) , we feel that we now find some evidence 
for further mass segregation in the difference between the 
LFs of unevolved stars in the Int and Ext regions of the 
cluster. 

4.2. Low mass star depletion in NGC 288 

In Fig. 6 the WF-Int and the WF-Ext LFs of NGC 288 
are compared with the LFs of M10, M22 and M55 by 
PZ99. The transformation from apparent to absolute mag- 
nitude has been obtained by assuming (to — M) = 14.67 
and E(B - V) = 0.03, according to Ferraro et al. (1999). 
The faintest point of the WF-Ext LF reaches M v ~ 10 
while the LF of the (more nearby) PZ99 clusters reaches 
My ~ 12-13. The faintest reported point appear to coin- 
cide with the turn-over point found in the LFs of the PZ99 
clusters, so this possible feature in the LF of NGC 288 is 
not detectable in our data (but it have been apparently 
detected by Pasquali, Brigas & De Marchi 2000, hereafter 
PBD00). 

On the other hand, the comparison in the range 6 < 
My < 10 shows that the LF of NGC 288 is significantly 
flatter than those of M10, M22 and M55. Note that the 
LFs of these clusters were obtained near the half-mass ra- 
dius (usually approximated with the half light radius r^, 
see Djorgovski 1993). This is the region where the mass 
distribution of the cluster population is expected to be 
fairly representative of the (present day) global popula- 
tion (see Meylan & Heggie 1997; Vesperini & Heggie 1997, 
and references therein). Since the Int and Ext sample are 
separated at 1 r^, the LF at the half-mass radius must be 
intermediate between the WF-Int and the WF-Ext, and 
the derived LFs should be representative of the global LF 
of NGC 288. 

The comparison presented in Fig. 6 strongly suggests 
that the global population of NGC 288 is significantly de- 
pleted of low mass stars with respect to M10, M22 and 



M55 (see PZ99 for an indirect comparison with other clus- 
ters). 

NGC 288 has a mass ranging from ~ i to <~ i that of 
the comparison clusters. Its central density ranges from 
~ \ t° ~ W (Pryor & Meylan 1993) that of the compar- 
ison clusters. Thus NGC 288 is intrinsically less resistent 
to harrassment and heating by external forces. Furthe- 
more, among the clusters considered here, it moves on the 
orbit having the smallest perigalactic distance, the high- 
est eccentricity and the highest inclination to the Galactic 
plane. Hence, it is expected to suffer strong bulge and 
disk shocks, and indeed may be one of the clusters with 
the highest disruption rates (Dinescu, Girard & van Al- 
tena 1999; Gnedin & Ostrikcr 1997). The coupled effects 
of internal energy equipartition (mass segregation), push- 
ing the lighter stars to outer, shallowly bound orbits, and 
of the tidal field of the Galaxy, removing the less bound 
stars venturing in the outer parts of the cluster halo, may 
well be at the origin of the depletion of lighter stars in 
NGC 288. We regard this result as a reassuring confir- 
mation of the theoretical predictions that account for the 
orbital characteristics of actual GGCs (Dinescu, Girard & 
van Altena 1999; Gnedin & Ostriker 1997). 

Recently, PBD00 presented a LF in the J and H near- 
infrared bands obtained with the HST NIC3-NICMOS 
camera (parallel observations). The sample is small but 
reachs fainter magnitudes (lighter masses) than our data. 
PBD00 conclude that the mass function of NGC 288 is 
very similar to that of all other studied clusters, and that 
it shows no sign of the strong harrassment predicted by Di- 
nescu, Girard & van Altena (1999) and Gnedin & Ostrikcr 
(1997). We regard the results (and the conclusions) of the 
present paper as more robust than those by PBD00 since: 
(a) the LF by PBD00 is based on a sample of 75 stars, 
a very small sample that may be subject to strong statis- 
tical fluctuations, and is not supported by artificial star 
experiments, while our LF is based on a sample of more 
than 6,000 stars corrected for all the observational biases 
with the artificial star experiments; and (b) the field ob- 
served by PBD00 is a small stamp (51'/2 x 51'.'2, <~ 15 times 
smaller than the field observed here) located at r ~ 2.4 77,, 
far outside the region of our sample which is expected to 
be representative of the global cluster LF. 

5. SUMMARY AND CONCLUSIONS 

The LF of the globular cluster NGC 288 (down to 
My ~ 10) has been obtained from a field covering a region 
comprised between the center of the cluster and r~2r(,. 

A new method to correct the observed LFs for all the ob- 
servational effects have been introduced and applied. The 
method is based on the Equivalent Sample concept, which 
may have many interesting applications in the study of 
stellar populations. 

The effect of the presence of binary systems in the final 
LF has been quantified and it has been found negligible, 
in the considered magnitude range. 

Comparison of the LFs obtained in different regions of 
the clusters, indicates the presence of mass segregation 
within NGC 288. Independently of the assumed normal- 
ization the LFs of the inner (WF-Int) and outer (WF-Ext) 
sample turns out to be different at the 99.1 % confidence 
level. 
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Fig. 6.— Comparison between the WF-Int and Wf-Ext LFs of NGC 288 with the LFs of M10, M22 and M55 by PZ99. The encircled points 
are the first points with Cf < 0.50. The horizontal bars indicate the normalization interval. 
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The comparison with the LFs of other clusters strongly 
suggests that the global population of unevolved stars 
in NGC 288 is significantly depleted of low mass stars 
(M < 0.5-0.6 Mq), in general agreement with the predic- 
tions of recent theoretical predictions taking into account 
also the orbital properties of the cluster. 
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Table 1 



PC-Int and WF-Int Luminosity Functions 







PC-Int 






WF-Int 






V 


N ohs 


Cf 


N c 




^obs 


Cf 


N c 






18.25 


15 


0.978 


14.8 


3.9 


126 


0.956 


135. 


.1 


12.2 


18.75 


29 


0.963 


30.2 


5.7 


264 


0.981 


270. 


.2 


17.1 


19.25 


44 


0.957 


46.3 


7.1 


376 


0.975 


388. 


.6 


20.8 


19.75 


70 


0.946 


73.0 


8.9 


433 


0.972 


452. 


.1 


22.8 


20.25 


56 


0.920 


60.7 


8.3 


519 


0.969 


544. 


.5 


25.1 


20.75 


51 


0.898 


58.5 


8.5 


522 


0.956 


551 


.5 


26.2 


21.25 


54 


0.832 


65.1 


9.1 


528 


0.947 


568. 


.6 


26.9 


21.75 


52 


0.767 


67.8 


9.8 


471 


0.935 


508. 


.3 


26.7 


22.25 


44 


0.695 


63.1 


10.0 


465 


0.927 


504. 


.9 


27.1 


22.75 


41 


0.607 


66.2 


10.9 


439 


0.913 


479. 


.9 


26.3 


23.25 


35 


0.490 


73.4 


13.3 


460 


0.899 


521. 


.2 


28.8 


23.75 


31 


0.203 






464 


0.796 


539 


.6 


27.5 


24.25 


2 


0.016 






322 


0.410 


813. 


.1 


52.3 


24.75 





0.000 






75 


0.096 








25.25 





0.000 






8 


0.011 









Note. — The V columns reports the V m (Vtruc) magnitude of the center 
of the 0.5 mag bins. N ^ s is the observed number of stars in the bin. Cf 
is the completeness factor in the center of the bin. N c is the corrected 
number of stars in the bin (i.e. this column contains the final corrected 
LF) and ctnc is the error on N c . The final corrected LFs are reported only 
down to the first bin with Cf < 50%. 
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Table 2 

PC-Ext and WF-Ext Luminosity Functions 





PC-Ext 


WF-Ext 




V 


N Q bs Cf N c a Nc 


^V obs Cf N c 





18. 


.25 


4 





.996 


18. 


.75 


13 





.989 


19. 


.25 


10 





.987 


19. 


.75 


18 





.980 


20. 


.25 


16 





.980 


20 


.75 


21 





.980 


21 


.25 


22 





.964 


21 


.75 


24 





.947 


22 


.25 


13 





.915 


22. 


.75 


15 





.895 


23. 


.25 


16 





.879 


23. 


.75 


23 





.806 


24. 


.25 


25 





.369 


24. 


.75 


5 





.069 


25. 


.25 








.004 



3.9 


2.0 


19 


12.7 


3.6 


70 


10.5 


3.4 


96 


18.5 


4.4 


116 


16.4 


4.2 


123 


21.8 


4.9 


162 


22.4 


5.0 


153 


25.3 


5.3 


144 


14.2 


4.1 


134 


16.3 


4.3 


132 


19.9 


5.2 


180 


25.6 


5.8 


177 


63.5 


14.0 


263 



204 
51 






.945 


19.4 


4.5 





.991 


70.8 


8.6 





.984 


97.4 


10.2 





.982 


118.6 


11.4 





.983 


124.2 


11.6 





.975 


168.9 


14.0 





.971 


157.8 


13.1 





.961 


149.2 


12.8 





.954 


143.1 


12.9 





.953 


137.3 


12.7 





.943 


192.2 


15.6 





.938 


173.4 


15.5 





.842 


299.7 


20.5 





.466 


419.9 


32.2 





.116 







Note. — The V columns reports the V m (Vtruo) magnitude of the center 
of the 0.5 mag bins. N \^ a are the observed number of stars in the bin. 
Cf is the completeness factor in the center of the bin. N c is the corrected 
number of stars in the bin (i.e. this column contains the final corrected 
LF) and ctnc is the error on N c . The final corrected LFs are reported only 
down to the first bin with Cf < 50%. 



